Black hole mass decreasing due to phantom energy accretion 
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Solution for a stationary spherically symmetric accretion of the relativistic perfect fluid with an 
equation of state p{p) onto the Schwarzschild black hole is presented. This solution is a generalization 
of Michel solution and applicable to the problem of dark energy accretion. It is shown that accretion 
of phantom energy is accompanied with the gradual decrease of the black hole mass. Masses of all 
black holes tend to zero in the phantom energy universe approaching to the Big Rip. 
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Observations of distant supernovae and cosmic mi- 
crowave background anisotropy indicate in favor of the 
accelerating expansion of the Universe . In the frame- 
work of a general relativity it means that a considerable 
part of the Universe consists of the dark energy: the com- 
ponent with a positive energy density p > and with a 
negative pressure p < — (l/3)p. This dark energy may be 
in the form of a vacuum energy (cosmological constant 
A) with p — —p, or a dynamically evolving scalar field 
with a negative pressure (quintessence or fc-essence 
13 )• One of the peculiar feature of the cosmological dark 
energy is a possibility of the Big Rip Q : the infinite ex- 
pansion of the universe during a finite time. The Big Rip 
scenario is realized if a dark energy is in the form of the 
phantom energy with p + p < Q. In this case the cos- 
mological phantom energy density grows at large times 
and disrupts finally all bounded objects up to subnuclear 
scale. 

What would be the fate of black holes in the universe 
filled with the phantom energy? Below we find the solu- 
tion for a stationary accretion of the relativistic perfect 
fiuid with an arbitrary equation of state p{p) onto the 
Schwarzschild black hole. Using this solution we show 
that the black hole mass diminishes by accretion of the 
phantom energy. Masses of all black holes gradually tend 
to zero in the phantom energy universe approaching to 
the Big Rip. The diminishing of a black hole mass is 
caused by the violation of the energy domination condi- 
tion p + p > Q which is a principal assumption of the 
classical black hole 'non-diminishing' theorems The 
another consequence of the existence of a phantom en- 
ergy is a possibility of traversable wormholes 6] . 

We model the dark energy in the black hole gravi- 
tational field by the test perfect fluid with a negative 
pressure, an arbitrary equation of state p{p) and energy- 
momentum tensor T^i, = [p + p)u^Ui, — pg^iu, where p 
and p are the dark energy density and pressure corre- 
spondingly, and = dx'^ / ds is a fiuid four- velocity with 
u^u^ = 1. The known analytic relativistic accretion so- 
lution onto the Schwarzschild black hole by Michel 
(see also 0) is not applied for a general case of a dark 



energy accretion. We adjust the Michel solution for the 
problem of dark energy accretion by elimination from the 
equations the particle number density. 

The integration of the time component of the energy- 
momentum conservation law T^^'j^ = gives the first in- 
tegral of motion for the stationary spherically symmetric 
accretion (the relativistic Bernoulli or energy conserva- 
tion equation): 
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where x = r/M, u — dr/ds, r is a, radial Schwarzschild 
coordinate, AI is a black hole mass and Ci is a constant 
determined below. To obtain an another integral of mo- 
tion we use the projection of the energy-momentum con- 
servation law on the four- velocity u^T^^.j, = 0, which for 
a perfect fiuid is 
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The integration of ^ gives the second integral of motion 
(the energy flux equation): 
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where u < in the case of inflow motion and a dimen- 
sionless constant A > 0. Note that the second integral 
of motion is obtained without use of the particle con- 
servation law. From ^ and one can easily obtain: 
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where C2 = ~Ci/A = Poo+p{poo)- From and Q one 
can find the relations for the fiuid velocity uh = u{2M) 
and density pn = p(2M) at the black hole horizon r = 
2M: 
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The constant A which determines the accretion flux is 
calculated by fixing parameters at a critical point. This 
provides us with the continuity of the solution from the 
infinity to a horizon. Following Michel we fix param- 
eters of critical point a; = a;, : 
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where c'^{p) = dp/ dp is the square of sound speed. Using 
(0) and l@J one can find the following relation: 
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which determines the density at a critical point p■^, = 
pixsf). Then for a given p*, with the help of lO, one can 
find a;* and u*. Constant A is fixed by substituting of 
the calculated values in Q . Note that there is no critical 
point outside the black hole horizon (a;* > 1) for < 
or Cg > 1. This means that for unstable perfect fiuid 
with < or > 1 a dark energy fiux onto the black 
hole depends on the initial conditions. This result has 
a simple physical interpretation: the accreting fiuid has 
the critical point if its velocity increases from subsonic to 
trans-sonic values. In a fluid with a negative or with 
> 1 the fiuid velocity never crosses such a point. It 
should be stressed, however, that fluids with c? < are 
hydrodynamically unstable (see discussion in 

Equations (O and along with the equation of state 
p = p{p) describe the requested accretion flow onto the 
black hole. These equations are valid for perfect fiuid 
with an arbitrary equation of state p = p{p), in par- 
ticular, for a gas with zero-rest-mass particles (thermal 
radiation), for a gas with nonzero-rest-mass particles and 
for a dark energy. For a nonzero-rest-mass gas the couple 
of equations Q and Q is reduced to similar ones found 
by Michel ^ . 

The black hole mass changes at a rate M — —Airr'^TJ 
due to the fluid accretion. With the help of (O and Q 
this can be expressed as 
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From this equation it is clear that the accretion of a phan- 
tom energy with poo +p(,Poo) < is always accompanied 
with the diminishing of the black hole mass. This result 
is valid for any equation of state p — p{p) with p + p < 0. 
If we neglect the cosmological evolution of poo then from 
© we obtain: 



As a particular fully solvable example we consider the 
perfect fiuid with a linear equation of state: 



p — a{p — po), a = const, po = const, 
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where a = c^ and it is supposed that p > 0. Among 
others this model includes radiation {p ~ p/3), ultra-hard 
equation of state [p — p) and the simplest models of dark 
energy (po = 0, a < 0). The constant a is connected with 
an often used parameter w — pj p by the relation w — 
a{p — po)/p- The physically reasonable case corresponds 
to < a < 1. 

In the linear model ^ the radius of critical point a;* = 
(l-|-3a)/2a, velocity at critical point ul = a/ (1-1- 3q;), ve- 
locity at the horizon mh = — (A/4)~"/'^^~") and constant 
A determining the energy flux in © is 
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This relation is valid only for stable fiuid with < a < 1. 
For unstable fluid with a < the constant A is indeter- 
minable by the condition of the solution continuity from 
critical point consideration. We suppose that in the un- 
stable case the growth of instabilities in the accretion flow 
will cause the asymptotic growth of the accretion velocity 
up to the limiting speed of light at the horizon uh — > — 1- 
This helps to fix the value of constant A: A = A. 

For some particular choices of parameter a the values 
p(a:) and u{x) can be calculated analytically. For exam- 
ple, for a = 1/3 the fiuid density is given by: 
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where 
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The density distribution for another physically interest- 
ing case a = 1 is given by: 
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where Mi is an initial mass of the black hole and a char- 
acteristic evolution time r = {AirAMilpoo + p{Poo)]}~^ ■ 



(11) 



The corresponding radial fluid velocity u = u{x) can 
be calculated by substituting of H10|) or (|ll|l into 
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For pq = the solutions H10() and describe corre- 
spondingly a thermal radiation and a fluid with ultra- 
hard equation of state. In the case of Poc < apo/(l + oi) 
the solutions I|1U|I and describe the phantom energy 
falling onto the black hole. For example, a phantom en- 
ergy flow with parameters a — 1 and pa = (7/3)/9oo 
results in a black hole mass diminishing with the rate 
M = -(87r/3)(2M)2p^. 

Now we turn to the problem of the black hole evolution 
in the universe with the Big Rip when a scale factor a(t) 
diverges at finite time (^] . For simplicity we will take into 
account only dark energy and will disregard all others 
forms of energy. The Big Rip solution is realized for 
p + p < and a < — 1. From the Friedman equations 
for the linear equation of state model one can obtain: 
|/?-fj3| oc a~^^^^°'\ Taking for simplicity po — we find 
the evolution of the density of a phantom energy in the 
universe: 

Poo = PocA ^1 - , (12) 

where 

.-'^-2afi(^....)"' (13, 

and Poo,i is the initial density of the cosmological phan- 
tom energy and the initial moment of time is chosen so 
that the 'doomsday' comes at time t. From © using 
(|12|l we find the black hole mass evolution: 

M = Mi{l + , 14) 

\ MoT T-tJ ' 

where Mq = (3/2) A^-'^jl + a\ and Mi is the initial mass 
of the black hole. For a — ~2 and typical value of A = 4 
(corresponding to uh = — 1) we have Mq = — 3/8. In 
the limit t ^ t (i.e. near the Big Rip) the dependence 
of black hole mass on t becomes linear, M ~ Mq (t — 
t). While t approaches to r the rate of black hole mass 
decrease does not depend on both an initial black hole 
mass and the density of the phantom energy: M ~ —Ma- 
in other words masses of all black holes in the universe 
tend to be equal near the Big Rip. This means that 
the phantom energy accretion prevails over the Hawking 
radiation until the mass of black hole is the Planck mass 

In remaining let us confront our results with the cal- 
culations of (not phantom) scalar field accretion onto the 
black hole 0, 0, Q, 0| . The dark energy is usually 
modelled by a scalar field </> with potential V{(f>). The 
perfect fluid approach is more rough because for given 
'perfect fluid variables' p and p one can not restore the 
'scalar field variables' 4> and S/(f>. In spite of the pointed 
difference between a scalar field and a perfect fluid we 



show below that our results are in a very good agree- 
ment with the corresponding calculations of a scalar field 
accretion onto the black hole. 

The Lagrangian of a scalar field is L = K — V, where 
if is a kinetic term of a scalar field 4> and F is a potential. 
For the standard choice of a kinetic term K = (j>-^^(j)'^ /2 
the energy flux is Tg^ = <j)^t4',r- Jacobson found the 
scalar field solution in Schwarzschild metric for the case 
of zero potential V ^ Q: cf) ^ ^^[t + 2Af ln(l - 2M/r)], 
where <f>oo is the value of the scalar field at the infinity. 
In |lj| it was shown that this solution remains valid also 
for a rather general form of runaway potential 1^('^). For 
this solution we have = — (2M)^(/)^/r^ and corre- 
spondingly M = 47r(2M)2(/)^. 

The energy-momentum tensor constructed from Jacob- 
son solution completely coincides with one for perfect 
fluid in the case of ultra-hard equation of state y = p 
under the replacement Poo ^lo/^i Poo ~* <^oo/2 UM- 
is not surprising because the theory of a scalar field with 
zero potential V{4>) is identical to perfect fluid consider- 
ation . In a view of this coincidence it is easily to see 
the agreement of our result (0) for M in the case oip = p 
and the corresponding result of [l^ . 

To describe the phantom energy the Lagrangian of a 
scalar field must have a negative kinetic term 0], for 
example, K = —(j>.^fj_(j>'^ /2 (for the more general case 
of the negative kinetic term see ,17]). In this case the 
phantom energy flux onto black hole has the opposite 
sign, Tor = —(t>^t(t>,r, where (f> is the solution of the same 
Klein-Gordon equation as in the case of standard scalar 
field, however with the replacement V — > —V. For 
zero potential this solution coincides with that obtained 
by Jacobson for a scalar field with the positive ki- 
netic term. Lagrangian with negative kinetic term and 
V{4>) = does not describe, however, the phantom en- 
ergy. At the same time, the solution for scalar field with 
potential V{(j)) = is the same as with a positive con- 
stant potential Vq = const, which can be chosen so that 
p — — (/)-^/2 + Vb > 0. In this case the scalar field repre- 
sents the required accreting phantom energy p > and 
p < —p and provides the decrease of black hole mass with 
the rate M = -A-K{2Mf4>l^. 

The simple example of phantom cosmology (without a 
Big Rip) is realized for a scalar field with the potential 
V = m^(\p- 12, where m ^ 10~^^ eV [ll|. After short tran- 
sition phase this cosmological model tends to the asymp- 
totic state with H ~ m(/)/3i/2 ^nd ~ 2m/3i/2. In 
the Klein-Gordon equation the term (with the men- 
tioned replacement V — s- — V) is comparable to other 
terms only at the cosmological horizon distance. This 
means that the Jacobson solution is valid for this case 
also. Calculating the corresponding energy flux one can 
easily obtain M = -47r(2M)20^ = -64M2to2/3. For 
Afo — Mq and m = 10^'^"^ eV the effective time of black 
hole mass decrease is r = (3/64)M~^77i^^ ^ 10"^^ yr. 

The possible physical interpretation of a black hole 
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mass diminishing is that accreting particles of phantom 
scalar field have a total negative energy [19|. The simi- 
lar negative energy particles are created in the Hawking 
radiation process and participate the Penrose black hole 
rotation energy extraction mechanism. Formally saying 
the black hole mass decreasing in the process of the phan- 
tom energy accretion is due to the violation of the energy 
domination condition in phantom energy. It should be 
noted that the existence of the horizon is not crucial for 
the decrease of the black hole mass due to phantom en- 
ergy accretion. The one-way-membrane property of the 
horizon merely makes this accretion inevitable and irre- 
versible. Any object or device would diminish its mass if 
it is capable to capture phantom energy. 
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